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Abstract 

We investigate the late time acceleration with a Chaplygin type of gas 
in spherically symmetric inhomogeneous model. At the early phase we get 
Einstien-deSitter type of solution generalised to inhomogeneous spacetime. 
But at late stage of the evolution our solutions admit the accelerating nature of 
the universe. For a large scale factor our model behaves like a ACDM model. 
We calculate the deceleration parameter for this anisotropic model, which, 
unlike its homogeneous counterpart, shows that the flip is not synchronous 
occurring early at the outer shells. This is in line with other physical processes 
in any inhomogeneous models. Depending upon initial conditions our solution 
also gives bouncing universe. In the absence of inhomogeneity our solution 
reduces to wellknown solutions in homogeneous case. We have also calculated 
the effective deceleration parameter in terms of Hubble parameter. The whole 
situation is later discussed with the help of wellknown Raychaudhury equation 
and the results are compared with the previous case. This work is an extension 
of our recent communication where an attempt was made to see if the presence 
of extra dimensions and / or inhomogeneity can trigger an inflation in a matter 
dominated Lemaitre Tolman Bondi model. 
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Introduction 



Following the high redshift supernovae data in the last decade [1] we know that 
when interpreted within the framework of the standard FRW type of universe (ho- 
mogeneous and isotropic) we are left with the only alternative that the universe is 
now going through an accelerated expansion with baryonic matter contributing only 
five percent of the total budget. Later data from CMBR studies [2] further corrob- 
orate this conclusion which has led a vast chunk of cosmology community ([3] and 
references therein) to embark on a quest to explain the cause of the acceleration. 
The teething problem now confronting researchers is the identification of the mech- 
anism that triggered the late inflation. Workers in this field are broadly divided 
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into two groups - either modification of the original general theory of relativity or 
introduction of any mysterious fluid in the form of an evolving cosmological constant 
or a quintessential type of scalar field. But as discussed extensively in the literature 
(we are sparing the readers here to repeat once again those arguments) both the 
alternatives face serious theoretical problems. In this context one important thing 
should not escape our attention. One intriguing fact in the framework of the stan- 
dard FRW model is that the accelerating phase coincides with the period in which 
inhomogeneities in the matter distribution at length scales < 10 Mpc become signif- 
icant so that the Universe can no longer be approximated as homogeneous at these 
scales. One should also note that homogeneity and isotropy of the geometry are not 
essential ingredients to establish a number of relevant results in relativistic cosmol- 
ogy. One need not be too sacrosanct about these concepts so as to sacrifice basic 
physics (energy conditions, for example) in relativistic cosmology. Conversely, if the 
universe is not apriori assumed to be homogeneous and isotropic, the observational 
data do not necessarily imply an accelerating expansion of the universe, or even if 
the cosmic expansion is accelerating it does not necessarily point to an existence of 
a dark energy. Thus to account for the observational data without introducing the 
concept of dark energy, varied arguments regarding the effects of inhomogeneities 
have been made and naturally a vast community of cosmologists have embarked 
upon a sort of 'mission' to explain (sometimes with conflicting claims) the observa- 
tional findings within inhomogeneous models. The immediate generalisation of FRW 
spacetime is the wellknown LTB model [1] which is also spherically symmetric but 
the spacetime is inhomogeneous. However, the assumption of spherical symmetry 
requires a centre of the universe so that the observer be located not too far from the 
centre to avoid undetected large anisotropy (a detailed study of LTB and allied cos- 
mologies and its relevance to current astrophysical issues may be found in [5]). The 
sojourn to the inhomogeneous path has a chequered history. Naively speaking there 
are two such arguments. One is that the apparent acceleration of the cosmology 
can be regarded as a result of an almost spherically symmetric but inhomogeneous 
peculiar velocity field, assuming that we are located at the vicinity of the symmetry 
centre [6, 7]. With this argument the acceleration of the cosmic volume expansion 
is not necessary. The other argument is that the acceleration of the universe is a 
physical reality and results from the backreaction effects due to the inhomogeneities 
in the background FRW universe [8, 9]. This idea is later supplemented by Carter 
et al [10] where the observed universe is assumed to be an underdense bubble in an 
Einstein-de Sitter universe and it was shown that from observational point of view 
their results become very similar to the predictions of ACDM model. However, 
in a recent communication Bolejko and Andersson [11] have calculated the volume 
deceleration based on Buchert averaging scheme and back reaction in some LTB 
models and have shown that for realistic cases the deceleration parameter turns out 
to be positive. At this stage a very brief mention of the formalism may not be 
out of place. The difference between the evolution of homogeneous models and an 
inhomogeneous universe is caused by backreaction effects, due to non linearity of 
Einstein equations such that the solutions for a homogeneous matter distribution 
leads in principle to a different description of the universe than an average of an 
inhomogeneous solution to the exact Einstein equations. So either we have to fall 
upon on exact solutions or invoke averaging the backreaction terms. If simply a vol- 



ume averaging is considered then such an attempt leads to Buchert equation. The 
Buchert equations are very similar to Friedmann equations except for the backreac- 
tion term, which is, in general nonvanishing if inhomogeneity is present (for a lucid 
review of the averaging scheme the reader is referred to [12, 13]). Moreover, the 
validity of the perturbative ansatz is questionable in that the claimed acceleration 
is later shown to be due to the result of extrapolation of a specific solution to a 
regime where both the perturbative expansion breaks down and the constraints are 
violated [1 1]. On the otherhand Kai et al [15] showed that if one hypothesizes the 
coexistence of expanding and contracting regions in space, the speed of the cosmic 
volume expansion can be accelerated. These models are constructed by replacing 
the spherical domains from the Einstein-deSitter model with a LTB dust sphere 
having the same gravitational mass. Another interesting suggestion has recently 
come from the works of Wiltshire et al [16] where a timescape cosmology has been 
proposed as a viable alternative to homogeneous cosmologies with dark energy. It 
realises cosmic acceleration as an apparent effect that arises in calibrating average 
cosmological parameters in the presence of spatial curvature and gravitational en- 
ergy gradients that grow large with the growth of inhomogeneities at late epochs. 
The model is based on an exact solution to a Buchert average of the Einstein equa- 
tions with backreaticon. Some people attempted to look into the problem from a 
purely geometric point of view - an approach more in line with Einstein's spirits. 
For example, Panigrahi et al [17] have recently toyed with the idea of dimension 
driven acceleration in a number of publications, where the extra terms coming from 
the higher dimensions create a sort of back reaction to drive inflation. Good thing 
about it is that one need not have to invent any exotic, unphysical matter field in 
this case. In an interesting contribution Wanas [18] introduced torsion to explain 
late acceleration. It is shown that torsion generates a new type of energy to be called 
torsion energy which is repulsive in nature, thus mimicing a quintessential type of 
field. While torsion inspired inflation has several desirable features(for example, 
geometrical origin) the problem with Wanas' model is that the geometry no longer 
remains Riemannian. 

Again we know [19] that in a matter dominated nonrotating model where parti- 
cles interacting with one another move along geodesic lines it is always possible to 
define a coordinate system which is at once synchronous ((700 — 1) and co moving. 
With this input Hirata and Seljak [20] claimed to have proved from Ray Chaudhuri 
equation [21] that in a perfect fluid cosmological model that is geodesic, rotation- 
free and obeys the strong energy condition (p + 3p) > 0, a certain generalisation 
of the deceleration parameter, g 4 must be always non-negative. But even with the 
perturbation considered by Kolb et al [8] the vorticity vanishes and consequently 
Kolb's claim is flawed. On the other hand, Iguchi et al [6] did obtain simulated 
acceleration in Lemaitre Tolman (LT) models with A = that obey the conditions 
set by Hirata et al, which subsequently led Vanderveld et al [7] to draw attention 
to this apparent contradiction between these two conclusions and to suggest that 
LT models that simulate accelerated expansion also contain a weak singularity, and 
in that case the derivation of HS breaks down. In addition to this, there are other 
singularities that tend to arise in LT models, and Vanderveld et al [7] have failed to 
find any singularity-free models that agree with observations. However in a pioneer- 
ing contribution Krasinski et al [22] neatly summed up the apparent controversies 



as also the claims and counter claims of different workers in this field and showed 
that the so called weak singularity is not a singularity at all, while the other types 
of singularities like shell crossing or shell focussing, generic to all inhomogeneous 
collapse may be taken care of with suitably chosen arbitrary functions appearing 
in the theory. Moreover, one should point out at this stage that unlike the ho- 
mogeneous case it is always difficult, if not a little ambiguous to define uniquely a 
deceleration parameter for inhomogeneous models. In this context Krasinski et al 
[22] also showed that Hirata-Seljak's formulation of is wrong, based on inadmis- 
sible averaging of the redshifts over directions and with the averaging dropped one 
gets correct signature of the deceleration parameter which may be both positive or 
negative. On the other hand Hansson et al [23] argued that when taking the real, 
inhomogeneous and anisotropic matter distribution in the semi-local universe into 
account, there may be no need to postulate an accelerating expansion of the uni- 
verse at all despite recent type la supernova data. In fact inhomogeneous structure 
formation may alleviate need for accelerating universe. 

In the present work our goal is completely different. In a recent communication 
one of us [2 1] examined the possibility if the presence of inhomogeneity or extra 
dimensions separately or jointly in an LTB model can achieve late acceleration 
without the aid of any extraneous scalar field. We found that while dimensions 
have no perceptible effect on nature of evolution the radial or angular acceleration is 
possible even in pure dust distribution if any of them decelerates fast enough in LTB 
model. The present work is an extension over that in the sense that we here introduce 
a Chaplygin gas as input. The lack of information regarding the provenance of 
dark matter and dark energy allows for speculation with economical and aesthetic 
idea that a single component acted in fact as both dark matter and dark energy. 
The unification of these two components has risen considerable theoretical interest, 
because on the one hand the model building has become considerably simpler and 
on the other hand such unification implies the existence of an era during which the 
energy densities of dark matter and dark energy are remarkably similar. Moreover, 
at present it is unclear whether the backreaction effects of inhomogeneities can 
actually accelerate the cosmic volume expansion. 

One possible way to achieve that unification is through a particular k-essence 
fluid, the Chaplygin gas with the exotic equation of state. While literature abounds 
with work on Chaplygin gas in FRW models ([25] and references therein), barring a 
few [26] we are not aware of works of similar kind directed to inhomogeneous space- 
time. However, relevant to mention that Gorini et al [27] have recently discussed 
a Tolman-Oppenheimer-Volkoff type of bounded distribution in presence of a gen- 
eralised Chaplygin gas. Here we discuss the evolution of a spherically symmetric 
inhomogeneous model with a Chaplygin type of matter field and get the interest- 
ing result that an initially decelerating phase transits to a late accelerating one in 
line with the current observational results. We compared our findings with those 
obtained via RayChaudhury equation also and get identical results. Organisation 
of our paper is as follows. In section 2 we solve the field equations for our inho- 
mogeneous spacetime with Chaplygin gas as matter field and have addressed the 
problems for both early and late time inflation assuming a flat 3 space under differ- 
ent subtitles. While our solutions are amenable to both early deceleration and late 



acceleration we get an interesting result that under suitable initial conditions the 
model also admits a bouncing type of universe avoiding the big crunch. In line with 
inhomogeneous collapse the bounce occurs at different instants for different shells 
unlike the synchronous homogeneous case. In section 3 we compared our findings 
with the conclusions coming also from the RayChaudhury equation. The paper ends 
with some concluding remarks in section 4. 

2 Field Equation 

ds 2 = dt 2 - X 2 dr 2 - i? 2 (r, t) dX 2 (1) 
where dX 2 represents a 2-sphere with 

dX 2 = d6 2 + sin 2 6d<f) 2 (2) 

and the scale factor, R(r, t) depends both on space and radial coordinates (r, t) 
respectively. A prime overhead denotes d/dr and a dot denotes d/dt. As pointed out 
in the last section the gauge, goo = 1 follows for a spherically symmetric, irrotational 
inhomogeneous system only when, p = 0. But in our case, p = only under extremal 
condition and not all through its evolution. So the ansatz, g 00 = 1 should be treated 
as an additional assumption to simplify field equations [28]. 

A comoving coordinate system is taken such that u° = l,u l = (i = 0, 1,2,3) 
and g^ u u^u v = 1 where Ui is the 4- velocity. The energy momentum tensor for a 
dust distribution in the above defined coordinates is given by 

^ = (P + P)«-P^ (3) 

where p(r, t) is the matter density and p(r, t) is the pressure. The fluid consists of 
successive shells marked by r, whose local density is time-dependent. The function 
R(r, t) describes the location of the shell marked by r at the time t. Through an 
appropriate rescaling it can be chosen to satisfy the gauge 

R(0,r) = r (4) 

The independent field equations for the metric (1) and the energy momentum 
tensor (3) are given by 

_ 2XR 1 + R 2 J_ f2R" R' 2 _ JCFT] _ 
G ° " XR + R? X 2 \ R + R 2 1 XR J ~ 9 (5j 

Gl= 4 +1 ~if-^ = - p (6) 
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Solving Gqi = equation we get 



XM) = if) (9) 
where f(r) is an arbitrary function of r. 

Since the WMAP data [13] shows that the universe is spatially flat to within a 
few percent we can take / = 1 such that the field equations finally reduce to the 
following two independent equations as 

R 2 R! R , 
R 2+2 R>R= P < 10 > 

2- + ^ = -p (11) 
R R 2 K 1 

From the the Bianchi identity we get for the inhomogeneous model the conser- 
vation law 

V V T^ = (12) 

which, in turn, yields 

S.p + -^S v [V=g(p + P )uV] + rO"u A = (13) 

V-g 

For Tq = and y/—g = X(r,t)R 2 (r,t)sin8, we obtain 

At this stage we assume that we are dealing with a Chaplygin type of gas obeying 
an equation of state 

A , *. 

P= 15 

P 

where A is a positive constant. It was first introduced as a cosmological fluid 
unifying dark matter and dark energy by Kamenshchik et al [25] and since has 
been widely studied in this context. Moreover it has found applications in particle 
physics via string theory [29] and its supersymmetric extension [30]. With this input 
we finally get 

'+^s(™ 2 )(H) =0 (16) 

which integrates to 
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Plugging in the expression of p from equation (10) we finally get 
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Figure 1: T/ie variation of R(r, t) and t for different values of C{r) is shown. The graphs 
clearly show that acceleration increases for greater C(r) i.e., greater inhomogeneity. 



This is the key equation in most accelerating models dealing with a Chaplygin gas 
in homogeneous models except that C (r) is now not a true constant but depends on 
space for inhomogeneous expansion. Moreover the above expression is not amenable 
to any closed form analytic solution but integration results in a hypergeometric 
series. Figure-1 shows that acceleration depends on C(r), which represents the 
inhomogeneity. So C(r) may be the measure of inhomogeneity in our case. 

CASE A: At the early stage of the cosmological evolution when the scale factor 
R(r, t) is relatively small the second term of the last equation (19) dominates and we 
get a sort of dust dominated universe for a particular value of C(r) = ^jV 4 yielding 

R(r,t)=r[t + to(r)]* (20) 

It has not also escaped our notice that with this scale factor we get a vanishing 
pressure when used in equation (11). Moreover for isotropic expansion (X' = R) we 



get p ~ as in FRW universe. Relevant to point out that the expression (20) is 
not exactly Tolman-Bondi like and our line element reduces to 



ds 2 = dt 2 -r 2 [t + t (r)]3 (dr 2 + rW) (21) 

If, at this stage, we assume that t (r) vanishes or losing its space dependence becomes 
a true constant (in that time translation is required) then we get 



ds 2 



dt 2 -r 2 tz(dr 2 + r 2 d{l 2 



(22) 



This is a new solution and may be termed as the generalised Einstein- deSitter 
metric for the inhomogeneous spacetime. In the analogous homogeneous case for 
zero pressure dust with vanishing spatial curvature(fc = 0) we get for FRW metric 
the wellknown Einstein-deSitter metric as(i? ~ is). From equation (18) we get the 
expression of density as 
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CASE B : (R(r, t) is very large) 
Type - 1: 

In the late stage of evolution the second term of the RHS of the equation (19) 
vanishes and we get 
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This yields a solution 



R(r, t) = Rq exp 



(24) 



(25) 



This is the wellknown de Sitter type of solution generalised to inhomogeneous 
space. A% behaves as A, the cosmological constant. However, if we make a radial 
coordinate transformation 



Ro exp 
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(26) 



the metric reduces to 



d s 2 = dt 2 - exp ( 2 \l t j df 2 - f 2 exp ^2 y ^ t ] dtt 2 (27) 

Type - 2: The volume expansion rate for our metric is defined through the 
4-velocity of the fluid, u a as 



3H = u a a = u a , b g ab = u a , b h 



ab 



(28) 



where 



h ab = g ab + u a u b ( 2 g) 

As commented earlier in the discussion while the definition works perfectly well 
for a FRW like homogeneous distribution of matter it is always a bit ambiguous to 
define the deceleration parameter of an inhomogeneous anisotropic model because 
the relation (28) does not take into account the directional preference of the matric. 
For example, Tolman-Bondi has a preferred direction, being the radial one. We 
can still give an operational definition to the average volume acceleration of our 
model. For inhomogeneous model the directional preference need to be emphasized 
in the expression for expansion. We define a projection tensor t ab that projects every 
quantity perpendicularly to the preferred spacelike direction s a (and of course the 
timelike vector field, u a ) such that 

t ah = g ab + u a u b - s a s b = h ab - s a s b (30) 

For our metric (1), 

f - ^/Sv (31) 

and the tensor projects every physical quantity in a direction _L to s a . One can 
now define the invariant expansion rates as 

H r = u a , b s a s b = ^ (32) 
H ± = -u a , b t ab = | (33) 



n " ,u R 



so that 



H = -H± + -H r (34) 
3 3 V ' 

Evidently the above definition gives a sort of averaging over the various directions 

for our anisotropic model. 

If one relaxes the condition of any particular preferred direction (like the radial one 
as in LTB model) one can explore the definition of the Hubble parameter in a more 
transparent way considering its directional dependence [31] as follows: 

H=\u a a + a ah J a f (35) 
3 

where a a b is the shear tensor and J° a unit vector pointing in the direction of 
observation. For an observer located away from the centre of the configuration it 
gives for our LTB case 

f-fH < 36 > 

where 6 is the angle between the radial direction through the observer and the 
direction of observation. Naturally when the two directions coincide, 6 = we get 




H = H r and for 9 = | it is H = He. A definition of deceleration parameter in a 
preferred direction can also be given in terms of the expansion of the Luminosity 
distance Dl in powers of redshift of the incoming photons. For small z one gets 

« = -* f S t + 1 < 37 > 

For = and 9 = ir /2 the acceleration is respectively 
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(38) 
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We shall subsequently see in section 3 that deceleration parameter defined this 
way has an important difference from what we later get in equation (56). Here the 
parameters do not depend solely on local quantities as opposed to the acceleration 
parameter of (56). For example we get via field equation(lO-ll) 

Mir) 1 , . 

where M(r) is the mass of the fluid distribution upto the comoving radial coordinater. 
Thus the equation (40) tells us that here the deceleration parameter q± depends on 
the total mass function and not on the local energy density of (48). 

One can look into the above expression of q± from a different standpoint also to 
assume a particular form of deceleration parameter as 



1 R a-R" 
q± = -HlR = bTR^ (41) 



where a, b and m are constants. Straight forward integration of equation (41) yields 

R(r, t) = R smh n u{r + t) (42) 



where, n = — , a = (Ro) n (- — l) and b = (Ro) n such that we get from equation 
(41) 

1 — n cosh 2 u;(r + t) 
n cosh 2 w(r + 1) 

showing that the exponent n determines the evolution of q±. A little inspection 
shows that (i) a < 0, i.e., n > 1 gives acceleration, (ii) a > 0, i.e., < n < 1 gives 
the desirable feature of flip, although it is not obvious from our analysis at what 
value of redshift this flip occurs. 

For n = 3, equation (24) is satisfied for A = ^ou 4 and in this case 

3 

q± = -sech 2 w(r + t) — l (44) 



Figure -2 shows the variation of q± and t for different values of r. We have seen from 
the graph that flip (i c ) occurs early at greater value of r, i.e., acceleration depends 
on inhomogeneity. The flip time (t c ) will be in this case 
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Figure 2: The variation of q± and t for different values of r is shown. The graphs clearly 
show that flip (t c ) occurs early at greater value of r, i.e., acceleration increases for greater 
r i.e., greater inhomogeneity. 



Type-3: 

It also follows from equation (18) that for the late universe (R ~ oo) 

< 4e) 



This is a mixture of a Cosmological Constant y A with a type of matter obeying 
a 'stiff fluid' equation of state. However it should be pointed out that it is an 
inhomogeneous and anisotropic generalization of the well known FLRW situation 
characterized by X(r,t) = R(r,t) where the quantities depend on time only 

Again as R(r,t) — > oo, we asymptotically get p = — p from this Chaplygin type 

of gas, which corresponds to an empty universe with Cosmological constant y^f • 
CASE C : 

Now we are trying to solve the equation (19) using the method of separation of 
variables. Let R(r,t) = g{r)a{t). From equation (19) we get 



where B = ir^— Constant (say). 

The equation (48) gives the hypergeometric solutions of a(t) with t. The solution 
and other features are same like homogenous case [25, 32] such as 

i) When A = 0, we get pressureless equation of state. Our solution reduces to 

2 

FRW type and in this case a(t) ~ frs. 

ii) At early stage of evolution, i.e., for small value of a(t), the equation (48) 
reduces to 3^- = ^ and we get a{t) ~ t§. 

iii) At the late stage of evolution, i.e., a(t) is large in this case, the equation 
(48) becomes (neglecting higher order terms ) 

3 ^ = VA + ^=a^ (49) 
a 2 y^A 

Solving the equation (49) we get the solution, 
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Figure 3: The variation of a(t) vs t is shown in this figure. Taking A=l & B=12. Flip 
occurs at t = t c = 0.661768 
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The nature of R(r, t) with t for a typical r is shown in the fig. - 3. 

In this case, R{r,t) = at t = ^A~^ In (j^j) ■ Here R = at t 

t c , which is the filp time. 
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CASE D : (For Negative Integration Constant) 

Considering negative value of C{r) presents interesting possibilities 
in that case the energy density increases with scale factor mimicing the phantom 
dark energy model and finally ending up as a cosmological constant. We get from 
equation(18)that for well behaved matter field the condition R' 2 R A > ^-jp, i.e. need 

to be satisfied. So a minimal value of scale factor exists for a typical value r, 

i 

which is R(r,t) min = 6 , pointing to a bouncing universe at early times. We 

thus see that the Chaplygin gas model interpolates between dust at small R and 
a cosmological constant at large R, but choosing a negative value for C(r), this 
quartessence idea lose. Following Barrow [34] if we reformulate the dynamics with 
a scalar filed ( and a potential V to mimic the Chaplygin cosmology, we see that a 
negative value for B dictates that we transform ( = i^. In this case the expressions 
for the energy density and the pressure corresponding to the scalar field show that 
it represents a a phantom field. This implies that one can generate phantom-like 
equation of state from an interacting generalized Chaplygin gas dark energy model 
in LTB universe. This feature has been discussed in the past [12]in the context of an 
effective description of inhomogeneous model evolving like a homogeneous solution 
following an averaging technique and also in details in [35, 27] . In our case the 
bounce is inhomogeneous in the sense that each shell characterised by a constant 
radial coordinater bounces at its own time. So the bounce is not synchronous each 
shell sharing a local dynamics. 



3 Raychaudhuri Equation : 

It may not be out of place to address the situation discussed in the last section with 
the help of the well known Ray Chaudhuri equation [21], which in general holds for 
any cosmological solution based on Einstein's gravitational field equations. The Ray 
Chaudhuri equation is 

0,X = % " 2 ( ff2 " ^) " l ^ + RuaV U V a (53) 

where the terms have their usual significance. With matter field expressed in terms 
of mass density and pressure Ray Chaudhury equation is finally given by, 

= -2(a 2 - c 2 ) - V - ^ (p + 3 P ) (54) 

in a co moving reference frame. Here p is the isotropic pressure. 

With the help of equation(34) we get an expression for effective deceleration 
parameter as 

H + H 2 
q = = -1 - 3 p (55) 

which allows us to write, 

6 2 q = 6a 2 + l2nG(p + 3p) (56) 
With the help of the equations (15), (18) & (56) we finally get, 



9 2 q = 6a 2 + 12nG 
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In our case the shear scalar evolves as 
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CASE A : Early Stage: At the early phase of this evolution when the scale 
factor R(r, t) is small enough the above equation reduces to 



2 q = Qa 2 + 12-kG 



R'R 2 



(59) 



It follows from the equation (59) that q, the deceleration factor is always positive. 
So accelerated expansion is absent in this dust dominated phase though inhomogene- 
ity is present here. Interestingly this result is very similar to the work of Alnes et 
al [36]. 

CASE B : Late Stage : 

Type - I: If we consider the late stage of evolution i.e., R(r,t) is large enough 
in this phase, the second term of the RHS of the equation (19) vanishes and we get 
from equation (57), 



9 2 q = 6a 2 - 24ttGVA 



(60) 



At this stage if we consider the scale factor given by equation (42) {n = |) the 
shear scalar becomes a 2 = |w 2 cosech 2 [2tu (r + 1)\ & A = ^§-oj 4l . The equation (57) 
reduces to 



9 2 q = 16w 2 cosech 2 [2cu (r + t)} - 32ttGuj' 



(61) 
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Figure 4: The variation of a 2 vs t is shown in this figure. Taking uj = 1 & r = 1. 



In figure - 4 shows a 2 vs t for a particular r. In this graph we have seen that 
as t increases a 2 decreases, i.e., when t — > oo, a 2 — > 0. So initially we get the 
decelerating universe and after flip it becomes accelerating in line with current 
observational result (see equation (60)). 

Type - II: Again if we consider first order approximation of equation (57), 
neglecting higher order terms, we get 



Q 2 q = 6a 2 + 2AnGA 



-A + 



C(r) 
R' 2 R 4 



(62) 



Let R(r,t) = g(r)a(t), so in this case a = which follows from the equation (58) 
Now the equation (62) reduces to 



6 2 q = 2AttGA 



a 6 



(63) 
Now q < 0, 



It follows from the equation (63) that flip occurs when a(t) = (-^) 

at a(t) > (-|) 6 i.e., acceleration takes place in this case. 

Negative Constant : If C(r) < 0, then B < 0, the equation (63) then becomes 



6 2 q = -2AttGA 



A 



B 



(64) 



From the above equation we have seen that always q < 0, which means we get always 
accelerating universe. 



4 Concluding Remarks 

The present work may be looked upon as an extension of one of our recent pub- 
lications where we examined the possibility in a higher dimensional LTB model if 
the inclusion of extra space jointly with inhomogeneity can induce late inflation in a 
dust model. While total volume acceleration is ruled out we found that preferential 
acceleration in radial direction is possible if the angular direction decelerates fast 
enough or vice versa. Here we have taken a Chaplygin type of gas as matter field to 
work out the same problem in a 4D spacetime. While there is a proliferation of work 
in the literature on homogeneous FRW model with Chaplygin gas we have not much 
come across work of similar type in inhomogeneous spacetime. Given the fact that 
it is difficult, if not a little confusing to define uniquely a deceleration parameter in 
inhomogeneous, anisotropic model we have nevertheless got the following definitive 
results. 

1. Aside from space dependence the mathematical structure and its followup in 
the section 2 is essentially similar to the works of homogeneous spacetime except 
the appearance of the term, C(r) in equation(17), which unlike its homogeneous 
counterpart is not a true constant but depends on the space coordinate. Its pres- 
ence introduces all the differences in cosmic evolution. Like FRW models our field 
equations are amenable to exact solution only at extreme values. We find that at 



early stage our solution reduces to inhomogeneous analogue of the Einstein-deSitter 
type of solution. 

2. In line with current observational findings our model accounts for early de- 
celeration and late acceleration showing the desirable phenomenon of flipin all the 
examples we examined. But as expected here the time of flip depends on space 
coordinate also. So flip is not synchronous as in FRW cases, occuring at different 
shells at different instants. So flip here is local, not global. Moreover flip occurs 
early for larger V. So for our spherically symmetric model the outer shells will 
start acceleration earlier and this is also a good news vis a vis when posited against 
the problem of shell crossing singularity generically associated with inhomogeneous 
models. 

3. Another interesting situation discussed is the possibility of bounce of our 
model from a minimum when the arbitrary function of integration, C(R) assumes 
a negative value. The bounce also shares the inhomogeneous characteristic of our 
model, the different shells characterised by r-constant hypersurfaces bounce at dif- 
ferent instants. Moreover we have here taken the original Chaplygin gas in our 
analysis but now generalised type of equation of state [37] are being increasingly 
used with greater freedom. In our future work we try to extend this work with these 
modified Chaplygin gas equations. 

4. To end the section a final remark may be reemphasised regarding the apparent 
accelerated expansion of the universe. To explain the SNIa observations the concept 
of accelerated expansion of the universe need to be invoked only for a FRW type of 
model. But one should point out that the Luminosity distance- Redshift relation, 
not the accelerated expansion is the quantity that can be directly measured. And 
within inhomogeneous models one gets better fit without the need to introduce the 
local accelerated expansion and consequent hypothesis of any extraneous, unphysi- 
cal matter field with large negative pressure. 
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